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This is a summary of a detailed study of heavy-light meson excited state energies. Our lattice 
measurements include both radial and orbital excitations. Particular attention is paid to the spin- 
orbit splittings, to see which one of the states (for a given angular momentum L) has the lower 
energy. In nature the closest equivalent of this heavy-light system is the B s meson, which allows 
us to compare our lattice calculations to experimental results (where available) or give a prediction 
where the excited states, particularly P-wave states, should lie. 
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1. Introduction 

The excited state spectrum of D s and B s mesons has attracted a lot of interest lately. For 
example, last year (2006) CDF and D0 collaborations reported measurements of two B s meson 
P-wave states [jjj], and BaBar and Belle measured excited D s meson states [^, Lattice QCD has 
now an excellent opportunity to offer some knowledge in the matter from the theory side. We have 
done an in-depth study of a heavy-light meson excited state energy spectrum on a lattice. This 
proceedings paper offers a summary of our results. Full details of the study can be found in m\. 

We have measured the energies of both angular and first radial excitations of heavy-light 
mesons. Since the heavy quark is static, its spin does not play a role in the measurements. We 
may thus label the states as L± = L ± i, where L is the orbital angular momentum and ±i is the 
spin of the light quark. 

Our main measurements are done on a 16 3 x 32 lattice with 160 configurations. The two 
degenerate quark flavours have a mass that is close to the strange quark mass (about l.lm s ). The 
lattice configurations were generated by the UKQCD Collaboration using lattice action parameters 
j8 = 5.2, csw = 2.0171 and K = 0.1350. The lattice spacing is 0.11 fm. More details of the lattice 
configurations used in this study can be found in Refs. ^ ||]. Because the light quarks are heavier 
than true u and d quarks, the pion mass is m n = 0.73(2) GeV. Two different levels of fuzzing (2 and 
8 iterations of conventional fuzzing) are used in the spatial directions to permit a cleaner extraction 
of the excited states. 

We also introduce two types of smearing in the time direction. First we try APE type smearing, 
where the original links in the time direction are replaced by a sum over the six staples that extend 
one lattice spacing in the spatial directions ("sum6" for short). To smear the static quark even more 
we then use hypercubic blocking ("hyp" for short), but again only in the time direction. The label 
"static" is used to denote the heavy quark that is not smeared in the time direction. 

2. Energy spectrum 

To obtain the energy spectrum we measure the 2-point correlation function 

C 2 (T) = (P t TG q (x,t + T,t)P !+T T' i U Q (x,t,t + T)) , (2.1) 

where U@(x,t,t + T) is the heavy (infinite mass)-quark propagator and G q (x, t + T,t) the light anti- 
quark propagator. P t is a linear combination of products of gauge links at time t along paths P and 
r defines the spin structure of the operator. The (...) denotes the average over the whole lattice. 
A detailed discussion of lattice operators for orbitally excited mesons can be found in In this 
study, the same operators are used as in [gj. The energies (mi) and amplitudes (a,) are extracted by 
fitting the C2 with a sum of exponentials, 

[Ci{T)]f\ j 2 « £ a Lfl e- miT a i)h , where AUx = 2-4, T < 14. (2.2) 

i=l 

In most of the cases 3 exponentials are used to try to ensure the first radially excited states are 
not polluted by higher states. Also 2 and 4 exponential fits were used to cross-check the results 
wherever possible. Indices f\ and f% denote the amount of fuzzing. 
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Figure 1: Energy spectrum of the heavy-light meson. Here L+(— ) means that the light quark spin couples 
to angular momentum L giving the total j = L ± 1 /2. The 2S is the first radially excited L= state. The 
D+— is a mixture of the D— and D+ states, and likewise for the FH — . Energies are given with respect to 
the S-wave ground state (IS). Here a = 0.110(6) fm is used to convert the energies to physical units. The 
error bars shown here contain statistical errors only. 
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2.805(14) 
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4.82(10) 
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2P- 


7.31(16) 


6.45(6) 


6.11(15) 


1P+ 


5.23(10) 


4.12(8) 


3.63(10) 


2P+ 


7.83(13) 


6.62(8) 


6.09(16) 


1D± 




5.48(10) 


5.22(17) 


2D± 




6.88(4) 


6.94(14) 


1D- 


6.08(14) 


5.10(7) 


4.79(12) 


2D- 


8.44(10) 


7.50(4) 


7.24(11) 


1D+ 


6.15(20) 


5.22(9) 


5.06(7) 


2D+ 


8.38(15) 


7.59(5) 


7.41(4) 


1F± 


7.21(12) 


6.09(10) 


5.69(8) 


2F± 


9.16(3) 


8.18(4) 


7.88(2) 



Table 1: Heavy-light meson energies on the lattice in units of ro. Here ro/a = 4.754(40) (+2 — 90) (Ref. [g]), 
which means ro « 0.52 fm. The n denotes the radial excitation and n— 1 gives the number of nodes in the 
wavefunction of the state. 



The extracted energy spectrum is shown in Fig. |T[ In most cases, using different smearing 
for the heavy quark does not seem to change the energies significantly — the exceptions being 
the P+ (and excited DH — ) state. The energy of the DH — state had been expected to be near the 
spin average of the D— and D+ energies, but it turns out to be a poor estimate of this average. 
Therefore, it is not clear to what extent the FH — energy is near the spin average of the two F-wave 
states, as was originally hoped. 

2.1 Interpolation to the b-quark mass 

One way to obtain predictions of the B s meson excited state energies is to interpolate in 1 /ttiq, 
where niQ is the heavy quark mass, between the static heavy quark lattice calculations and D s 
meson experimental results, i.e. interpolate between the static quark (mg = oo) and the charm quark 
(niQ = m c ). Here we, of course, have to assume that the measured D s meson states are simple 
quark-anti-quark states. This is not necessarily true: for example the mass of the D* (2317) is 
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Figure 2: Interpolation to the b-quark mass. The ratio m c /mi, is taken to be 0.30(2) (from [g]; shown 
by the vertical band). The D s meson experimental results are from [|[ (blue circles), and the B s meson 
experimental results are from [||] (blue circles) and [Jl[] (green triangles). TrinLat group's m c /m.Q = results 
(maroon diamonds) are from [|To|]. Our results (using "hyp" configurations) are marked with red squares. 
The agreement with the measured B s meson P-wave states is very good as seen in Table ^. If the "sum6" 
lattice results were used instead of the "hyp" results, the predicted energies at m c /mh would be too high 
compared to experiments. 





hyp 


sum6 


experiments 


0+ 


405 ± 20 MeV 


399 ±12 MeV 




1+ 


447 ± 20 MeV 


441 ±12 MeV 




1+ 


477 ± 22 MeV 


5 17 ±22 MeV 


460 ±3 MeV 


2+ 


489 ± 22 MeV 


528 ± 22 MeV 


470 ±3 MeV 



Table 2: Our predictions for B s meson mass differences, M(B*)-M(B V ), for the P-wave states. The uncer- 
tainty in the ratio m c /nib was not taken into account in the error estimates. The experimental results are from 
[0, §]. Here the "hyp" lattice P+ energy is taken to be E(P— ) + spin-orbit splitting = 440(30) MeV instead of 
the direct lattice measurement E(P+) = 423(39) MeV. The latter would lead to predictions E(l + ) = 466(28) 
MeV and E(2+) = 477(28) MeV for the B s meson. 
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Figure 3: Interpolation to the b-quark mass: higher excited states. The lines illustrate what the interpolation 
would look like, if the D s meson states were D-wave states. The experimental results are from [|J ||]. 
Interpolating to m c /mt, predicts D-wave J p = 3~, J p = l~ at 938(23) and 826(42) MeV respectively. 



much lower than what is predicted by conventional potential models, and it has thus been proposed 
that it could be either a four quark state, a DK molecule or a D% atom. However, the inclusion of 
chiral radiative corrections could change the potential model predictions considerably [11]. Here 
we assume simple quark-anti-quark states and do a linear interpolation (see Fig. ||). The two lowest 
P-wave states seem to lie only a couple of MeV below the BK and B*K thresholds respectively. Our 
predictions are given in Table ^. 

As for other excited states, BaBar and Belle observed two new states, D* 7 (2860) and D* sJ (2700), 
in 2006 [§, ||]. The J p quantum numbers of the D* 7 (2860) can be + , 1 , 2 + , etc., so it could be a 
radial excitation of the D* (2317) or a J p = 3~ D-wave state. The first interpretation is rather pop- 
ular, but our lattice results favour the D-wave J p = 3~ assignment in agreement with Colangelo, 
De Fazio and Nicotri [|l2|]. Interpolation then predicts a D-wave J p = 3~ B s state at 938(23) MeV. 
The alternative interpretation as a 2P state, when compared with our lattice result, would lead to 
an unreasonably large 1/niQ dependence. In addition, the D* 7 (2700) could be a radially excited 
S-wave state or a D-wave J p = 1 state. If the latter identification is assumed, then a D-wave 
J p = 1 B s state at 826(42) MeV is expected (see Fig. Ejh. 
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Figure 4: The spin-orbit splittings of P-wave and D-wave states. The P-wave spin-orbit splitting is small 
(consistent with zero), whereas the D-wave spin-orbit splitting is much larger and positive. 



E(nL+) - nE(L-) 


Indirect 


Direct 


E(1P+) - E(IP-) 
E(1D+) - E(ID-) 
E(2P+) - E(2P-) 
E(2D+) - E(2D-) 


-5 ±45 MeV 
102 ±54 MeV 
-10±83MeV 
70 ±60 MeV 


11 ±13 MeV 
124 ±18 MeV 



Table 3: P-wave and D-wave spin-orbit splittings using "hyp" configurations. 



2.2 Spin-orbit splitting 

We pay particular attention to the spin-orbit splitting (SOS) of the P-wave states, i.e. the 
energy difference of the 1P+ and IP— states. We extract the SOS in two different ways: 



1. Indirectly by simply calculating the difference using the energies given by Eq. 22 , when the 
P+ and P— data are fitted separately. 

2. Combining the P± and P— data and fitting the ratio C2(P+)/C2(P— ), which enables us to 
go directly for the spin-orbit splitting, «iip + — mip_. 

The D-wave spin-orbit splitting is also extracted in a similar manner. The results of the fits are 
given in Table |3| and in Fig. ||. The two smearings, "sum6" and "hyp", give different spin-orbit 
splittings for the P-wave: 74(10) MeV and 11(13) MeV, respectively. This is unexpected, and is 
studied in detail in A similar difference is not seen in the D-wave, where direct analysis gives 
124(18) MeV ("hyp") and 91(19) MeV ("sum6"). 



3. Conclusions 

• With the "hyp" lattice, our predictions for the 1 + and 2 + P-wave state masses agree very well 
with the experimental results. We also predict that the masses of the two lower P-wave states 
(0 + and 1 + ) should lie only a few of MeV below the BK and B*K thresholds respectively. 
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Also with the "hyp" lattice, the P-wave spin-orbit splitting is small (essentially zero), but 
the D-wave spin-orbit splitting is clearly non-zero and positive. In contrast, another lattice 
group finds the D-wave spin-orbit splitting to be slightly negative (see [|l0|]), i.e. they seem 



to observe the famous inversion [|13|]. On the other hand, in Qll| ] Woo lee and Lee suggest 
that the absence of spin-orbit inversions can be explained by chiral radiative corrections in 
the potential model. 
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